Dufour and Soret effects in a magnetized and non-magnetized 

plasma 



L. S. Garcia-Colm 

Departamento de Fisica, Universidad Autonoma Metropolitana-Iztapalapa, 
Av. Purisima y Michoacdn S/N, Mexico D. F. 09340, 
Mexico. Also at El Colegio Nacional, Luis Gonzalez Obregon 23, 
Centra Histdrico, Mexico D. F. 06020, Mexico. 

A. L. Garcia-Perciante 
Depto. de Matemdticas Aplicadas y Sistemas, 
Universidad Autonoma Metropolitana-Cuajimalpa, 
Av. Pedro A. de las Santos No. 84, Mexico D. F, Mexico 

A. Sandoval- Villalbazo 

Departamento de Fisica y Matemdticas, Universidad Iberoamericana, 
Prolongacidn Paseo de la Reforma 880, Mexico D. F. 01210, Mexico. 

(Dated: February 5, 2008) 



1 



Abstract 

It is well known that magnetic fields affect heat conduction in a different way in the 
direction parallel and perpendicular to the field. In this paper a formal derivation of this 
phenomenon and analytical expressions for the transport coefficients based in the Boltz- 
mann equation are presented. Moreover, the Dufour effect or diffusion thermo effect is 
usually ignored in plasma transport theory. In this work, analytic expressions for the paral- 
lel and perpendicular thermal conductivities as well as the coefficients for both the thermal 
diffusion, or Soret, effect and the Dufour effect are formally derived. It is also shown how 
the heat conduction in the perpendicular direction decreases with increasing magnetic field 
and how in both directions the diffusion thermo effect is significant compared to thermal 
conduction. Other aspects of this work are also emphasized. 
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I. INTRODUCTION 



The study of transport processes in plasmas is about to reach its seventieth birth- 
day. Indeed, in 1936-1937 Landau published his pioneering paper by casting the 
Boltzmann collision integral in the form of a Fokker-Planck equation thus regarding 
that collisions in plasmas may be treated diffusively. Using very intuitive arguments 
he then arrived to what is today known as the collision integral in the logarithmic 

. . n 

approximation for a gas with Coulomb interactions [Ifl. This method has become 
the basis of practically all approaches to the study of such processes in many ap- 
plications. For instance, most of the work involving astrophvsical problems use this 
method as further developed by Spitzer and collaborators jj, 1^ over fifty years ago, 
also later by other authors and independently by Braginski [6|, |7|. The last 

effort along these lines is the apparently not fully recognized but magnificent work 
on transport theory in plasmas written by Balescu in 1988 [8|. In our opinion it is 
probably the best and most complete report written on the subject as of today. After 
Balescu's work, research on the kinetic theory of plasmas seems to have diverted in 
two directions. One is related to the kinetic problems which arise in fusion-related 
research. In tokamaks and other devices, the geometry and inhomogeneities present 
in the relatively large magnetic fields that prevail, affect the particle dynamics in a 
crucial way. This topic is referred to as neoclassical transport [9 1 and the literature 



here is abundant, a sample of which may be found in Refs. [lO|, lUl and the works 



cited therein. The second application is characteristic in astrophysical problems 
where the Spitzer-Braginski approach seems to dominate. Typical situations are re- 



lated to the problem of thermal conduction in tangled magnetic fields 

^ nn 



12|, thermal 



conduction in clusters of galaxies [13| and similar problems [lJ|[l5|. A few papers 
deal with traditional kinetic theory but they do not seem to proliferate. An example 
of a formal kinetic treatment of a two temperature gas can be found in Ref. [16| 
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and of a multispecies plasma in Ref. [l7|. However, a common denominator in most 
recent works in plasma kinetic theory is the use of the Landau-Fokker-Planck version 
of the kinetic equation and few have used the full Boltzmann equation to study this 



ISj it is clearly 



question. Indeed, in a recent authoritative monograph on the subject 
stated that Braginski's equations "are now being universally accepted and because 
of their completeness and correctness are considered standard". We shall come back 
to these statements later on. 

On the other hand, another approach to the study of transport phenomena in 
plasmas, using the full properties of the Boltzmann equation, was first attempted 
by Chapman and Cowling in 1939 in the first edition of their monumental work on 
kinetic theory [l^. It is very likely that since they had in mind mostly applications 
to the physics of metals, their method did not permeate to other branches of physics 
as smoothly as the Spitzer-Braginski ones did. In 1960 W. Marshall wrote a series of 
three reports entitled "The Kinetic Theory of an Ionized Gas" in which he extended 
the preliminary results of Chapman and Cowling to the calculation of the transport 
coefficients of such a system when no magnetic field is present and in the presence of 
a magnetic field. Apparently these papers were never published in a regular journal 
and have thus been ignored in most work done in this subject 2o| . 

Another question that has been completely ignored in the several approaches to 
the transport processes, with the sole exception of Balescu's work js], is the deriva- 
tion of the basic assumptions of Classical Irreversible Thermodynamics (CIT) from 
a kinetic model. This not just means discussions about the approach to equilib- 
rium of a charged system, the H-theorem, and the explicit calculation of the entropy 
production exhibiting consistency with the second law of thermodynamics. It also 
includes the derivation of the complete equations relating the appropriate fiuxes with 
their corresponding thermodynamic forces, the calculation of transport coefficients 
involved therein, specially those associated with the so-called cross-effects, and then 



showing explicitly that the matrix of such transport coefficients is a symmetric one, 
the well known Onsager reciprocity theorem. Curiously enough, even for an in- 
ert multicomponent mixture of dilute gases this latter property was obtained only 



recently [2l|. This complete program has been accomplished for the so-called "vecto 
rial" transport processes in a dilute plasma, heat conduction, diffusion and electrical 
conduction. However the complete work, which includes viscomagnetism, is simply 
too bulky to attempt a publication in a single regular article ^2^. 

In this paper we would like to discuss with certain detail the problems of heat 
conduction and diffusion using the full Boltzmann equation as a formalism which 
not only involves the full two body dynamics of the interacting species but it also 
allows one to fulfill the program described in the previous paragraph. One reason for 
doing so is that these processes have become rather important in many astrophysical 
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^plications, specially in .he enipnatie p.oblen. of the so-called coding flows 



The main result which we want to emphasize here is that besides the 
heat transfer rates which arise from simple heat conduction due to thermal gradients, 
results which are essentially in agreement with those of Spitzer and Braginski, the 
influence of diffusion (mass transport) is rather surprising. With or without the 
presence of a magnetic field, the diffusion thermo effect or Dufour effect, as well as 
the pressure thermal effect, seem to play a rather important role in enhancing heat 
conduction. To our knowledge this particular result has never been published. It is 
as striking as the thermal analog of the well known Hall effect, the so-called Righi- 
Leduc effect, which we have reported in a separate paper [25|. Although similar 
results are also present in the case of electrical conduction we shall not deal with 
them here to keep the size of this paper within a reasonable length. 

Yet, another reason for carrying out this work is to provide the kinetic foundations 



of magnetohydrodynamics within the framework of CIT which, as mentioned above, 



was partially done in Ref. [8|| and later on extended to the relativistic regime 



The structure of the paper is as follows. In Section II we briefly summarize the 
tenets of classical irreversible thermodynamics, Section III is devoted to the kinetic 
theory behind the calculation. In Section IV we compute and discuss the transport 
coefficients and Section V is left for some concluding remarks. 



II. THERMODYNAMIC BACKGROUND 



Classical irreversible thermodynamics [271, l28| establishes that for any physical 



system in a non-equilibrium state the conservation equations which govern the time 
evolution of the corresponding local state variables together with the assumption that 
the local entropy is a time independent functional of the subset of such variables, 
which are even under time reversal, lead to the well known entropy balance equation, 
namely 

^ + V.J. = .. (1) 

In Eq. dl]) p and s are the local density and entropy respectively, Jg is a vector 
determining the entropy flux through the system's boundaries and a the entropy 
production which may be expressed as a sum of products of the fluxes Jj produced 
by the thermodynamic forces Xj (both quantities being, in general, tensors) 

^T = J]j.oX,, (2) 

where "o" indicates a product of equal rank tensors. 

To supply the sufficient information so that the conservation equations are well 
determined one assumes that linear constitutive equations hold true 

j 

where the matrix L^j, the transport coefficients matrix, must be symmetrical L^j = 
Ljk, the well known Onsager's reciprocity theorem. 



For the particular situation to be considered in this paper, a 
an inert ionized gas, Eq. (JS]) would be normally written as follows 
denote the heat and mass fluxes then 



Dinary mixture of 



3- 



If Jq and Jjui 



Jg = -LnVT - L.^dij , (4) 

— # — * 

Ji = -L21VT - L22dij , (5) 

where dij = —dji is the well known diffusive force. The coefficient Lu is the thermal 
conductivity, L12 is related to the diffusion thermo effect, L22 is the generalized 
diffusion coefficient and L21 the thermal diffusion coefficient well known as the Soret 
effect. More details about them will be discussed later on. In the dij representation 



it has recently been shown that for an ordinary multicomponent ideal mixture 



21I 



712 — ^21 , 



(6) 



which is Onsager's theorem. This theorem, as pointed out in Ref. [21| does not hold 
for arbitrary choices of the thermodynamic forces but we shall not pursue this matter 
here. 

The main objective of this paper is to derive Eqs. ((4|) and ^ for the case of the 
two component mixture present in an ionized gas, namely electrons and ions, in the 
presence of a magnetic field with explicit expressions for the transport coefficients and 
to carefully analyze their significance. The proof of Onsager's reciprocity theorem 



follows at once from the results of Ref. 



2l| in the absence of a magnetic field. With 



a magnetic field the proof will be given elsewhere. 



III. KINETIC THEORY 



In order to carry out the above program we make use of the Boltzmann equation, 
implying that our system is a dilute mixture of electrically charged particles with 
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charges Cj, masses rrii (i = a, h) and assume that the ions are protons {Z = 1). If 
fi is the single particle distribution function of species i defined as usual in /i— space 
(f, Vi) then, 

i,j=a 

Here J (fifj) is the well known collision kernel Q] taking into account collisions 
between particles of different and same species. Also, to avoid effects arising from the 
long range Coulomb potentials, the integration limits are calculated with a cutoff at 
a distance A^, the Debye cut-off distance which is the well known value of the impact 
parameter allowed for the definition of a collision. 

The electromagnetic force is introduced as the Lorentz force where the magnetic 
field is weak enough such that the prevailing Larmor frequency, namely uji = CiB /rrii, 
is such that tUjT (where r is the characteristic time between collisions) is not to large. 
Thus, we are considering a regime in which collisions prevail over cyclotron motion. 
In this force, the electric and magnetic fields are the self-consistent fields derived 
from the Maxwell equations Q]. 

Once Eq. ([7|) is clearly defined one ought to proceed with the derivation of the 
conservation equations and the proof of the H-theorem. The former is a standard 



step widely discussed in the literature 



30 



3l|. The second one is a little bit more 



22| but leads, for a wide class of 



tricky due to the presence of the magnetic field 
external potentials, to an equilibrium state described by a Maxwellian distribution. 
Here, for simplicity, we are assuming that Ta = Ti, = T. Also it is quite clear that 
the solution to the homogeneous part of Eq. ([7]) is given by 

fj^ (f, t)=n (f, t) [-^^^jp^^j e-a^ , (8) 

where Ci = Vi — u{r,t), u{r,t) = {riaUa + UbUb) /n are the chaotic and barycentric 
velocity respectively, n = Ua + Ub is the number density and all these quantities 
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are locally defined. The local temperature T (r, t) is defined through the internal 
energy density e (r, t) by the standard ideal gas law. In Eq. ([8|) the electric field has 
been ignored for simplicity, as will be in the rest of this article but can be easily 
incorporated in the formalism through its potential as a conservative force. 



We thus proceed directly to the solution of Eq. (l2l) f o 
Enskog method widely discussed in the literature |19l . 



lo wing the Hilbert-Chapman- 



30 



3l| . The method is based 



on two assumptions, firstly that fi can be expanded in a power series in Knudsen's 
parameter e which is a measure of the magnitude of the macroscopic gradients within 
the extension of the physical system and secondly, that the time dependence of fi oc- 
curs only through the local state variables, n (r, t), u{r, t) and e (r, t), or equivalently 
T(f,t). Then 



fi (f, Vi, t) = ff^'^ {r,Vi\ n, u, T) 



(1) 



r, Vi \ n, u 



(9) 



where n, u and T are functions of r and t. For e = 0, /j = ff^^ which when substituted 
in the conservation equations leads to the Euler equations of magnetohydrodynamics. 
To first order in e, since the right side of Eq. ([7]) is already of order e we get, after 
substitution of Eq. ([9]) in Eq. ((71) and a lengthy but straightforward manipulation of 
the left side, the linearized Boltzmann equation, namely 

b 

■Ci = T^^e 



nii 0^ 
— QQ : Vu + 



rriicf 



2kT 



VT n r 
, 1 da 

2 T m ^ 



irti 



Ci X B 



(1) 



dvi 



+ C{p 



pkT ^ 

j=a 



a, b , 



where the superscript "0" over a tensor indicates its symmetric and traceless part. 



form is also well known 
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„ ,"ni\ niUi 



n 



are the linearized collision kernels whose explicit 
sol . l3l| and dij is the diffusive vector force given by 



n 



TTli 



rrii 



P 



P 



niTij 
pp 



jriiCj 



rrijei) u x B , 



(11) 



which satisfies the property dij = —dji. 

The solution to Eq. ( fTOl) follows from the standard technique in the theory of 
integral equations using the fact that its structure resembles that of a system of 
linear inhomogeneous integral equations. Neglecting the term ^CiCi : Vm, which 



involves second rank tensors and is important only when studying viscomagnetic 
effects 22|, one gets that (see Refs. 21\ and js^), 

y^f ) = A, ■ VT + ■ 4- , (12) 

where Aj and Dj are the most general vectors that can be constructed from the base 
vectors in our space Cj, q x 5 and (^cl x x B = B (^cj ■ B^ — B'^Cj, whence 

Ki = + Af xB + Af^B (^Ci ■ B^ , (13) 



and 



^(1)^^ + of V, xB + 5 ( Q • 5 I , (14) 



where Af\ ©f^ {k = 1, 2, 3) are functions of all scalars in this space, 
(n, T, clB', (q-5)'). 

The next step is to substitute Eqs. ( fT3l) and ( fT4l) into Eq. ( ITOl) and separate terms 
in VT and (ijj. The algebraic procedure is both lengthy and tricky. It is outlined in 
Appendix A. The final result is that one gets a set of two linear integral equations 
for each of the functions Aj and Dj namely 

/f (0 - I) "5 = {C + C m, + c,K^)} , (15) 



+ {C (qA) + C (qA + c, A)} , (16) 
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where TZi = A 



(1) 



B^Af^ for i = a,b. Also, in Eq. ([M 



(17) 
(18) 



where G^"* and G^'' are defined in Appendix A. 



An identical set of equations is obtained for D. except that the inhomogeneous 
term is fi^^ ^ck, the Ai's are substituted by Vi's and ^ by a function /C also defined in 
the Appendix A. Equations f fTSl ) and ( |T6l) and their analogs for the Vi functions are 
the furthest results that can be obtained without resorting to the explicit form of the 
interaction potential between the different species. The evaluation of the linearized 
collision terms involved in the homogeneous parts of these equations will lead to the 



standard collision integrals 



20 



29|. 



The solution to Eqs. (jTSi l and (1161 ) is somewhat lengthy. We fist notice that 
Eq. (fTSi l does not contain any explicit contribution of the magnetic field, the TZi 
function does not involve any property other than those already associated with the 
Aj's so that its solution can be immediately implied from the one arising for the 
same problem involving only an inert dilute mixture and its solution is available in 



any standard book on kinetic theory 



3, 



30, Isjj , although we will come back to it 



in Appendix B. On the other hand Eq. ( ITGll is different on account of the first two 
terms in the right side involving the magnetic field (iB). To accomplish its solution 
we follow the almost orthodox methods in kinetic theory, namely to expand the un- 
known functions A^-'^'' (and D,-*^-*) in terms of a complete orthonormal set of functions, 
the Sonine (Laguerre) polynomials. Thus we write that 



A}' 



A; = 1, 2, 3 



(19) 



m=0 



where the Sonine polynomials are defined elsewhere 



29 



3l| and have the property 
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that 



J^^ e-^'S^\x) S^\x) x-^-Hx = "^^^^6,, , (20) 
and the first few of them are Sn\x) = 1, Sn\x) = —x + n + 1 and 

S^:\x) = i (n + 1) (n + 2) - (n + 2) X + ^ . (21) 
In Eq. (fT9ll the coefficients a'p'^"^'' are functions of the scalars n, T, c^, i?^, and 



Ci ■ Bj . A similar expression holds for 
Clearly, for Eq. dill), 

oo 

A = A« + ^5Af ) = J2 ^T^Stii {cf) , (22) 

with 

^M^^(i)M^^^^{2)M_ ^23) 

Equations (fTOl l. (p2ll and ( l23l) considerably reduce both the structure of the inte- 
gral equation [Eq. (fT6i l] and its analog for the 3^^^ functions as well as the subsidiary 
conditions, Eqs. f IClll in Appendix C. In fact, introducing the dimensionless velocity 



^* = \/^^,, (24) 



noticing that f-^'' = rii (2^^)^^^ ^ ^iiid using Eq. (l20ll we get, after straightforward 
integration that 

X:n.(ap)+i?\(^)(°))=0, (25) 



J]n,ar=0, (26) 



this implying that all coefficients a] for n > are not restricted, and also that 



6 = E^*^4°'. (27) 
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and a similar expression for the function /C defined in the Appendix A. 

As we shall see below, the solution to the integral equations (fTSl) and (fT6l) is 
considerably simplified if we first compute explicitly the fluxes in terms of the forces, 
Eqs. ([4]) and ([5]), and obtain formal expressions for the transport coefficients. This 
will be the subject of the following section. 



IV. TRANSPORT COEFFICIENTS 



This section is devoted to the explicit calculation of both the mass and heat fiuxes 
using the results of the previous sections and some rather well known results of 
irreversible thermodynamics. Let us begin with the transport of mass. The diffusion 
fiux for species a is defined as 

Ja = mJ Caf^^if^^Uca. (28) 



Substituting Eqs. (fTTl) and ffT2]) into Eq. (1281) and noticing that all the resulting inte- 
grals contain the dyad CaCa^ we get with the use of Eqs. ( fT9l) and (l20l ) that 



/oo 
;(0) ^ af^-^S^^^CaCjca = ^kTa^^m , (29) 
m=0 

where I is the unit tensor and use has been made of the fact that the contribu- 
tion arising from the symmetric traceless part of CaCa vanishes. The same result is 
obtained for the ©j functions so that after collecting all terms we get that, 

fa = Uak {a«(°) VT + 42)© VT xB + af ^ [b ■ Vt) B 

+T [d^:^^'^l, + xB+[b. I,) B } , (30) 

and a similar expression for Jb since Ja + Jb = and dab = —dba- Equation (i30l l is a 
very eloquent result namely, knowledge of mass fiux requires not the full information 
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contained in the infinite series of Eq. (11911 but only of six coefficients, three for each 
of the thermodynamic forces, VT and dij. This fact will considerably simplify the 
solution of the integral equations derived in the previous sections. Notice also how 
the magnetic field infiuences the diffusive, thermal and mechanical effects which give 
rise to the mass fiux. This is readily appreciated examining Eq. (ITTll which even 
in the absence of an external conservative field has two contributions, one arising 
from the standard diffusion (the first two terms) and another one arising from the 
magnetic force term u x B. Moreover, the thermal diffusion effect (proportional to 
VT) is given by the first three terms in Eq. ( l30ll , each one representing a different 
contribution. Referring our vectors to a cartesian coordinate system, if we choose B 
along the 2-axis, B = Bk, then the third term is simply B'^^k which applies as well 
to the last term, so that Eq. fl30l) may be rewritten as 



Ja = n„fc{(a«(°)+i?2af(°))V||T + a«WVxT + ai2)(o)5V,T 



+ T 



1)(0) 



+ B'S^^('^) 4,1, + dW(°)4,xT + S^^^'^Bl,, I , (31) 



where ViiT = ^k and such term represents the mass fiux parallel to the magnetic 



field; V_lT = + is the thermal force giving rise to the fiow of mass in the plane 
perpendicular to the magnetic field and finally V^T = — jr-i is the thermal force 



dz ' 
dx ' dy - 

dx dy 

producing a mass fiow which is perpendicular to both the direction of the field and 
to VT. Similar features also hold for the diffusion force d^j which will be analyzed 
in detail later on. We must also mention that when 5 = 0, this last contribution 
vanishes and Eq. (l3Til reduces to 

Ja = riak a«(°) VT + riakT S^^^'^d^t , (32) 

which is of the standard form required by non-equilibrium thermodynamics Icf. 
Eq. (jn])]. Here Uaka^a^^^'^ is the so-called Soret or thermal diffusion coefficient 27, 28|. 



In the presence of a magnetic field there will be three contributions to the Soret effect 
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which are easily written down from Eq. (ISB and the results of Appendix B. Indeed, 
if we denote these coeficients by S, then we get that 

5„ = ^ X 2.9 



2 smK 

where ip is defined in Eq. (IB16I) . Similarly, using Eqs. flB3l) and (|B19l) one gets that 

3^="^ (56.3 - 662x2 - 2.25x^) , 
2Ai ^ ' smK 

where 

Ai (x) = 19 + 2078x2 + 2650x^ + 9x^ . (33) 

Here x = cUgT where the mean free time r, is obtained directly from the collision 
integrals (see Appendix B), and is given by 

^i, ' 

Finally, using the same equations 

8, = "^ (647x + 2.2x=^) 



2Ai ^ ' smK 

The values of these coefficients as functions of x are plotted on Fig. 1. 

Let us now turn our attention to the fiow of heat. According to the tenets of CIT 



the heat current is defined as 



27|, 



4 = '^-E^^— ' (35) 
^ ^ rrik 

where hk is the enthalpy of the fc*'^ component which is equal to |fcT for an ideal 
gas. Thus in kinetic theory, the analog of this equation is clearly given by 



4 = ^.-|^^( — + -1 ' (36) 
2 \mamb' 
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or using the standard definition of Jq as well as Eq. (i28l) we get that, 



^;=*^i:/(i|-|)5//°vs"''5 (37) 

i=a ^ ' 

We now perform the following steps: substitute Eqs. (I12H14I) and (fTOll in Eq. (l37ll and 
after examination of the six resulting terms realize that all of them have as common 
factor the integral 

tji'^-l) l^f^' t ^"-'S^dc. = -lkTj: ^.f <" , (38) 

i=a ^ ^ m=0 i 

and a similar result for the d[''^ functions. This leads to the result that 

= _^fc2y^^|ap)vT + aP)5x VT + aP) • Vt) 5 

i * 

+ T [df^^''>dat + df^^'^B X dab + rff^'^ [b ■ dab) b] } , (39) 

where care has to be taken in the second summation to account for the property 
that dab = —dba- Equation ( l39ll contains all possible contributions to the flow of heat 
including those arising from the magnetic fleld since m x i? is present in the diffusive 
force dab- Moreover, each term has the same structure of the two contributions 
appearing in Eqs. ([30]) and (l3B so that there will be three directions of heat flow, 
the one parallel to B and which is unaffected by the fleld, one perpendicular to B 
and a third perpendicular to both B and VT. These latter one is the analog of the 
Hall effect in electromagnetism and is known as the Righi-Leduc effect. This effect. 



discovered in 1887, was carefully studied by W. Voigt in 1903 



elsewhere 



25|. 



28|, 



32| and discussed 



Clearly, if i? = the ordinary (Fourier) thermal conductivity is simply given by 

i=a 
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whose calculation demands only knowledge of two coefficients af'^^^\ i = a, b, not 
restricted by subsidiary conditions, instead of the full knowledge of the infinite series 
involved in Eq. (fT9l) . In the diffusion thermo effect, usually referred to as the Dufour 
effect, a little care is required to examine the transport coefficient since even if 
i? = 0, dab contains two contributions, one arising from V (^) and another one, 
the so-called pressure diffusion coefficient proportional to Vp [c.f. Eq. (JTTl)]. The 
resulting equation is of the form 

4 = «:VT - Dudab , (41) 

where 

Du = l{kTfi:^d^''''^ (42) 

i=a * 

which is of the standard form required by non-equilibrium thermodynamics. 

The final step to be undertaken here is to calculate the twelve coefficients which 
appear in Eqs. (|3T1) and (JSQ]), respectively and that requires the explicit solution of 
the integral equations (fTSl ) and ( flGl ) and their analogs for the ©•'^^ functions under the 
assumption stated in Eq. (fTOi l together with the subsidiary conditions, Eqs. f l25ll and 



(1261 ) and their analogs for the df'^ coefficients. This involves, once more, a rather 
lengthy and strictly mathematical procedure including the detailed evaluation of 
fourteen collision integrals which for completeness reasons, is outlined in Appendix B. 
In what follows we shall simply quote the results there obtained. Let us concentrate 
on the first line of Eq. ( l39ll . set B = Bk along the 2;-axis and write for the conductive 
part of the heat fiux, 

'j'Y = -KiiViiT-Kj_V±T-KsVsT, (43) 



where each term can be easily identified form Eq. (I39|l . Thus, using Eqs. (1B16I) and 
(IB17I 1 we obtain, 

5k'^T J 
Kii = X 2.01nr — — . 44 

4 irie sKm 
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As expected, k\\ does not depend on B and has been obtained taking into account 
that for Z = 1, a hydrogen plasma, nip = 1836me and assuming further that = 
Ub = n/2, the plasma is fully ionized, k is Boltzmann's constant, rrie and e are the 
electron mass and charge, respectively, eo the electrical constant in Farads — and 
tfj is Coulomb's logarithmic function 



■i/) = In 



(45) 



obtained from the collision integrals. Debye's characteristic shielding length, A^, is 
defined as 

A.= y^^. (46) 

It is important to mention that Eq. ( l44ll is in agreement with the Spitzer-Braginski 
result to within 10% as the interested reader may easily verify. On the other hand 
the expressions for k± and Kg are strongly field dependent and according to the set 



of equations (|B19l) are given by 



5 k J- Tl/T~ / _ o /I \ J / V 

= — X (38.7 + 2270a;2 + 161x^) — — , 47 

4 nie Ai sKm 

5 k'^T nr , ■i\ J / x 

Ks = — X (206x + 2644x^) — — , 48 

4 me Ai ^ ' sKm ' ^ ' 

and Ai is defined in fl33l) . Figure 2 clearly exhibits the behavior of the three conduc- 
tivities for weak fields. A number of comments about these results is pertinent. In 
the first place kj,, the thermal conductivity for the heat fiow along the (x, y) plane, 
reduces to n\\ when 5 = as may be easily verified. Since in this case /t^ = we 
recover Eq. (14T1) . that is the ordinary Fourier's constitutive equation. These results 
for the parallel and perpendicular components of the heat fiux (for cjr ^ 1) agree 
with those in Ref. Secondly, the Righi-Leduc conductivity [Eq. (l48ll ] vanishes 
when i? = and for a whole range of values B ^ Q never exceeds the value of n\\ . 
This is in agreement with Balescu's results 
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On the other hand, it is also interesting to compare with the original results 
obtained by Marshall [20||. To do so we must clarify exactly and precisely where the 
differences are. We begin with the equation for the electrical current written at the 
bottom of page 22 in Marshall's report number 3, which leads to Eq. (3.61). This 
correct result is identical to ours. Next, he defines a new diffusive force 

D = ^ cTi (49) 

eim2 - e2mi 

where di is given in Eq. (3.14). So far, so good. Then he writes the electrical current 
as shown in Eq. (3.62) and concludes that the coefficients a/, an and am are the 
thermal diffusion coefficients. This is questionable. Firstly, the force D has three 
components, a coefficient times V(^), a coefficient times Vp and the coefficient of -E', 
see Eq. (3.63). The coefficient of V(^), a concentration gradient that produces an 
electrical current is called the "Dorn effect". The coefficients of E' are the electrical 
conductivities, when we have the case Vp = (the pressure contribution to j). The 
coefficients 0^, and 0^^^, a thermal gradient contributing to j, are the Thomson 
coefficients (thermoelectricity) none of which have to do with thermal diffusion. As 
we shall point out later, mass and electrical currents are easily related to each other 
[see Eq. (l58ll ] whence, except for a numerical factor, ^j^^, the coefficients of VT, 
0^-0^^^, are indeed the same as the three Soret coefficients which we have explicitly 
obtained in a different context (see Eq. (l3Tll and Fig. 1 ). On the other hand, except 
for a numerical scaling factor, cr^, a^^ and a^^^ (Eqs. 7.8-7.10) are proportional to 
the ordinary diffusion coefficients not discussed here. 

Next we consider the equation for the flow of heat. In the form written by Mar- 
shall in Eq. (7.16), the 6*^, 6^^ and 6^^^ coefficients, the thermal conductivities given 
in Eq. (7.19) are not precisely our Ky, and because our definitions of heat flow 
differ by the term ^kT{ni (ci) + n2 (C2)) which according to CIT should be included 
in [see Eqs. (l35l) and ( l36l ) ] whereas Marshall deals with it a la Chapman and 
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Cowling leading to complicated expressions for the thermal conductivities. Never- 
theless the 6''s and k's have been here compared as shown in Fig. 4. It is important 
to emphasize that the expression used by us as the definition of the heat flux allows a 
clear separation between thermodynamic fluxes and forces and has become standard 



23|. 



since the classical work of de Groot and Mazur 

Finally, the coefficients Ci~Ciii in Eq. (7.20) are related to the so-called Benedicks 
effect in the presence of a magnetic field (coefficients of E') and the diffusion thermo 
effect (coefficients of V (^) and Vp), the standard Dufour coefficients, which we shall 
label "all together" and call the -Dy, D± and Dg. These are absent in Marshall's work. 

Indeed, in an entirely analogous fashion the second line of Eq. (I39l ) can be written 



as 



/ ^ ^ 

yJqj = -Dl\dab\\ - D_i_dab± " Dgdabs , (50) 

where using Eqs. flB20l) leads to, 

5 (kT)'^ 

D\\ = — X 0.29nr (51) 

" 4 me 

^ 5 (fcTf^ ^ ^5 g _ gg 23,2 _ Q 22a;4) (52) 
4 me Ai ^ ^ ^ ^ 

^ 5 (^ nr ^ 

4 me Ai ^ ^ ^ ^ 



Once more when B = 0, D± = D\\ and we recover Eq. f l4T]l identifying 1^12 with 
The magnitude of these coefficients are shown in Fig. 3. 

Let us examine in detail the structure of the diffusive force dij when the calcula- 
tions are performed in the comoving system so that u = 0. For the fully ionized gas 
V— = and therefore, 

n \ p J p 2 T 



20 



using also the fact that rrij ^ mi. For this mentioned above, the heat flux 

is enhanced by the thermal pressure difl^usion since Eq. (l39|) now reads 

4 = -^er {k|V||T + fi;lV±T + KlVsT] , (55) 

where = /« + ^ for each term. Note that, even in the absence of a magnetic field, 
the parallel thermal conductivity is enhanced by a factor of only about 7%. 

The magnitudes of the effective and standard thermal conductivities are shown in 
Figs. 4 and 5. 

To finish this section we must also emphasize the fact that there is a magnetic 
contribution both to mass and energy transport that arises from the term 

dii = ^ (rriiej — rrijei) u x B , (56) 
pp 

which also contains the here neglected effects of an electric field which, if present, 
will appear in Eq. (l56|l as E' = E + u x B . This modification is rather important 
when studying the problem of electrical conduction but we shall leave the details for 
a future publication. 

V. CONCLUDING REMARKS 

In this paper we have used the standard Boltzmann equation to study a fully 
ionized dilute plasma in a magnetic field. The main results obtained are explicit 
formulas for the transport coefficients associated with mass and heat fiux as well as 
their cross effects. These latter ones have been completely ignored in the literature; 
here we stress their importance for different values of the magnetic field for a given 
particle density and temperature under the weak field approximation so that if r is 
the mean collision time, UeT, where Ue is the Larmor frequency for electrons, is not 
too large. The full results are shown in the previous section. These results confirm 
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Balescu's conjecture stating that the values of the transport coefficients are not 
very sensitive to the method used in obtaining them. Also, we obtain an enhanced 
heat conduction due to the Dufour cross-effect which leads to an effective thermal 
conductivity in all directions [Eq. (I55D ] even in the absence of a magnetic field. 

We could have easily added an external field E without any effort. This would 
lead to the study of electrical conductivity since in the expression for dij, Eq. (fTTl) . 
the last term contains the full Lorentz force E+ux B. Since the electrical conduction 
flow is defined as 

b 

/e = V— 4, (57) 



rrii 

i=a 



and Ja + Jb = 0, 



Je = eJa , (58) 

marub 



thus, using the results for the mass flow, Jg can be readily computed from Eq. (l58l) . 
We avoided this discussion because of the intrinsic interest of electrical conduction 
which deserves a separate treatment. This has been discussed in another paper [34 1. 



The restriction here to a fully ionized hydrogen plasma is merely didactic. Taking 
Z ^ \ and setting = ?C (r, t) n where X (r, t) < 1 is a trivial step, all the 
collision integrals in the Appendix B can be written in terms of Xif^-C). In such 
a case 7^ and besides the pressure diffusion coefficient there will be a true 
Fickian contribution. The reader may pursue this generalization with no problem 
whatsoever. 

Finally, a word is indispensable concerning the results here obtained and those 
obtained by Balescu 8| and Braginski using the Fokker-Planck (FP) equation. 
Balescu uses the Grad-like moment expansion method to compute the transport 
coefficients both for electrons and ions in a plasma for different values of Z. As we 
mentioned in the previous section the values for the thermal conductivities of the 
electrons is rather satisfactory. No comparison is possible with any cross coefficients 
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since neither he nor Braginski calculate them. However, a word of caution must 
)e here introduced. Moment methods, as has been clearly pointed out recently 
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36l |. have to be dealt with carefully. On the one hand they may be subject to 



inconsistencies and secondly, when computing a transport coefficient, say with N (>5) 
moments, it is not clear which are the relevant contributions arising from different 
orders in the gradients. Although this problem may arise in the Chapman-Enskog 
expansion it is automatically taken care off by Knudsen's parameter. Therefore 
the different contributions in successive order in the gradients arise, Navier-Stokes, 
Burnett (second order), super-Burnett (third order) and so on. A fair comparison of 
the two methods would in principle require the introduction of such a parameter in 
the results derived from any moment expansion to classify the different contributions 
in terms of the gradient order. Despite its drawbacks, this method is still used in 
plasma kinetic theory calculations (see for example Ref. 3^ ) 

The results obtained by Braginski are carefully examined by Balescu so it is not 
necessary to repeat them here. Let us only add that in his result the relationship of 
the FP equation with the tenets of GIT is not pursued. The fact that the character- 
istic vector space to which the solution of the Boltzmann equation pertains is absent 
prevents the full discussion of the various cross effects when B 0. The so-called 
frictional force introduced in the solution of the FP equation is artificial, the macro- 
scopic (barycentric) velocity u is not a thermodynamic force and according to Curie's 



principle 



27l | the only vectorial fiows which may couple together are the heat, mass 



and charge fluxes. Thus we are convinced that our results are finer in their content 
than those arising from the FP equation. We will not pursue this discussion further 
since it would lead us away from the objective here acomplished. 



23 



Appendix A: INTEGRAL EQUATIONS FOR AND % 



In this appendix we outline the procedure leading to Eqs. (fTSl) and (fT6|l starting 
from Eq. (fTOl) . Substitution of Eq. (fT2ll into Eq. (fTOll yields two equations, 



.(0) / ^jCg _ 5 



J a 



(Al) 



+ n 



(0) 



C AJ +C Aa + A 



an identical equation for species b and two analogous equations for and Db which 
differ from Eq. (lAip in their inhomogeneous term, ^Cafa^\ Thus, the procedure to 
be applied to Eq. (lAll) is identical in the three other cases. Substitution of Eq. (fT3ll 
into Eq. (lAll) splits each term in this equation into three terms one for each scalar 
i = 1, 2, 3. Define M as 

^ ^^-[A«c; + Af (4x5)+Af5(c„-5; 
5^e,| dc,/f A, (c, X 5^ 



function 

M 



-/W^fc-x^). 



_f(0) 
J a 



pA;T 



(A2) 



and examine term by term. The first one 



CaX B 



) ■ ^ (A«4) = (4 X 5) AO 



since the second term vanishes. For the second one. 



B] Al^) 



(A3) 



CaX B 



_d_ 

dCa 



d 



Ai^) (^ca X e)] = A(2) (^4 X 5j ■ — (4 X 5 



+ 4 X B] ■ 



Cn X B 



(2) 



9A^ 
54 



(A4) 
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we make use of the vector identity 



Vc-, (caXB]= (VCa) X 5 = I X 5 , 



and thus 



CaXBj- [Ca xBj = [Ca X B j ■ [I X B j = [Ca X B j ■ B - B'Ca 

where the subscript reminds us that all derivatives are taken with respect to Qj. The 
third term in the first line of Eq. flA2l) vanishes. 

After exchanging (^Cj x B^ ■ Ca = — (^Ca x B^ ■ cj in the second line of Eq. ( 1A2I1 
we are left with 



f ^c, + A}^^ [cjxB]+ Af'B [cj ■ B 



.(2) 



,(3), 



■iCnX B 



(A5) 



and the first term yields 

e, j dc,ff^P^c,c, = \Y.^,j dc,ffp^pc^ . (A6) 



The second one is (\Y.j(ij ! dcjff^P^fc]i^ X B. However, [(^ 

Ca X B^ X B and since the third term vanishes, B ■ (^Ca x B^ = 0, we get after the 
use of the identity 



dcj 



(0), 



(A7) 



for any arbitrary even function of Cj, R{cj), that 



M 



B^Ca-B{ca-B 



B^Ca-B[ca-B 



G 



(2) 
B 5 



(A8) 
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where 



G 



(1) 

B 



G 



(2) 
B 



2 ^ 

3 
3 



dcjf. 



(0) A (1) 



1 



(0) a(2) 
3Jj 



dcjf. 



52 
1 



(A9) 



(AlO) 



When Eqs. (|A8l) - (|Al0l) are substituted into Eq. (lAip we get a linear form in terms of 
the basic vectors q,, x _B and B ^c^ ■ B^ . Setting the coefficients equal to zero we 
obtain three equations, namely, 

e.B^ 



2kT 2 



m„ 



^B^Gl + G (A«4) + G (a«c. + 



4 ''^b y ' 
(All) 



m 
-G 



Af ) f 0, X B 



Al^) I c„ X 5 



(A12) 



(A13) 



= -^A(^)4 + ^li^c-Gg) + G (Af c- ) + C (Ai=^)c- + Af ^c,) . 

Multiplying the third equation by B"^ and adding it to the first one yields Eq. (fTSll 
in the text whereas multiplying the second one by iB and adding it to the first one 
yields Eq. f flGl) in the text. Identical results are obtained for species b and also for 
the functions i = a, b; j = 1, 2, 3. The only difference, as we have pointed out 
above is in the inhomogeneous term, and instead of the G function we obtain 



K 



B 



1,2. 



(A14) 



Appendix B: COLLISION INTEGRALS^ 



As mentioned in Sect. IV, the evaluation of the transport coefficients requires the 
knowledge of twelve coefficients of the series given in Eq. (fTOl) and its analog for the 
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(k) 

functions, namely 

«M = a^^^""^ + B'a^^ t = a,b, m = 0, 1 , (Bl) 

^(m) ^ ^ ^2^mrn) , ^ ^ ^ m = 0, 1 , (B2) 

the m superscript denoting the order of the polynomial in Eq. (fTOll and 

^M^^(i)M^^^^(2)M^ (B3) 

with the same values for i and m. Notice that af^^ = a\^^ = af^^'^^ ii B = and 
the same relation holds for the d^s. 

The determination of the twelve independent coefficients in Ea s. (|Bll) -( lB4ll is 
achieved using a variational method 



wards used by several authors 



m 



irst introduced by Hirschfelder [29| and after- 



20 



22l |. It is based on the linearity property of the 



linearized Boltzmann collision kernel. Let Gij = Gij {wi, Wj) and Hij = Hij {tVi, Wj) 
be any two functions of the dimensionless velocities Wi and Wj. Define 

[G,j, Hi,\^ = — [ ... [ {H[^ - Hij) ft^ffg^.cr (Q) dQdwidw, , (B5) 

where the subscript ij in the bracket denotes integration over the variables Wi and 
Wj and the differential cross section a {VL) dfl = bdb de for a collision between two 
particles of species i and j with an impact parameter b{< A^) and orbit inclination 
e. 

Using the symmetry properties of the kernel in Eq. ( 1B5I1 . the fact that [ ]-j is a 
linear operator and for any two arbitrary functions Ki = Ki (wi) and Lj = Lj (wj) 
the quantity {K, L}, 

{K, L} = J2 niUj [Ki + K,, U + 1,^^ , (B6) 
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is such that 



{K, L} = {L, K} , 
{K, L + M} = {K, L} + {K, M} 



(B7) 
(B8) 



and 



{K,K}>0, (B9) 

where for obvious reasons the equahty follows if and only if is a linear combination 
of the collisional invariants. Equation (IB9I1 is the basis that sustains the variational 
procedure. Let Ri and Rj be the sought solutions to Eg . fflSTl. Let ti (tVi) {i = a, h) 



be a solution to this equation. Then one can prove that 



20 



22 



29| 



b 

/ dwih (wi) : U (wi) = {U, U} < {Ri, Ri} 



(BIO) 



where Ii{wi), the inhomogeneous term, is known. Thus the trial function t must 
be chosen so to maximize the collision integral {Ri, Ri}. For the trial function we 
propose that 

M 

^ — - (Bll) 



p=0 

SO that M characterizes the order of the approximation. The left side of Eq. (IBlOll 
can be readily integrated after substitution of Eq. ( IBllI ) and use of Eq. (i20i ). This 
leads to the result that 



15kT f + ^4^) = {to, t,} . 



(B12) 



Ki = Li = af'ci + ap^S'g/gCj 



We now choose that, to a first approximation, m = 0, 1 so that 

i = a,b. (B13) 
It is readily shown that for a binary mixture Eq. (IB6P gives 

{K, L} = 2nl [K,, L,l^ + 2nl [K^, U],, + 2n,nb{[K,, L^l, 

+ [Ka, LbU + [Kt, L^U + [Kt, . (B14) 
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Equations (IB13I1 and (IB14I1 . after a lengthy but straightforward algebra, lead to 
a system of algebraic equations for Qf\ ai^^ and whose coefficients are 
fourteen collision integrals which must be explicitly evaluated. This is an exercise, 
■;edious but rather straightforward, which yields the following results when mf, ^ rria 
2J, 



bb 







[Wb, Wb]ab = ^iV , 



2^' 



Wb, Wb\Wf,-- 



[Wa, WbU 



- ab 



Wi, Wi I w. 



i = a, b , 



1 W^a - ^ j , Wb [wl - ^- 



13 



ab 



15 



-Lp. 



MiLp. 



ab 
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ab 



Wa, Wb W^ 



' /\/f3/2 

-M;^ if 



'Wa\W^~-],Wb 



ab 



ab 



-2^1 



2 5\ ^ / 2 5 
Wa\W^ - -j , Waiw^ - - 



Wb\w^--j,Wb\w^-- 



(B15) 



bb 
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where 

rria + rrib 

and 

_ 1 _ 

^ = ^ " 4 (2vr)3/^ y^eg (fcT)^/^ ' ^ 
where all symbols are defined in Eq. (I44l) and in Eq. (145|1 . The mean free time 
r is defined here as the inverse of the collision integral y?; we stress here that the 
logarithmic function is a consequence of the evaluation of the collision integrals. 
Using these collision integrals and Eq. (IB16I) the resulting algebraic system, too long 
to be given here, is maximized to satisfy Eq. (IBlOll making use of the subsidiary 



conditions Eqs. (l25ll and (l26l) . After the system is solved one finds that 

= 2.94r = a«W + E^af W , 

= l.96r = a«(i) + B^af^^'> , (B17) 

a^^ = 107r = 4^)^^) + E^af . 

Here we have introduced the restrictive (but completely unnecessary) condition that 
the plasma is fully ionized, = = f so by Eq. (l25ll = a^^\ Also rrib ^ 
(mp = 1836 me) so that Mi ~ ^. 

(k) 

An identical calculation for the functions leads to the results that 
40) = = l.2r = + E^rff , 

5(i) = 0.29r = d«« + i?243){i)^ (B18) 
5« = 1.5 X 10-V = + . 

The solution to Eq. (fTBl) is slightly more sophisticated due to the presence of the two 
field dependent terms in its right side. The method followed is apparently due to N. 
Davison {2^, 33 1. One proposes a trial function Q and constructs a functional D (Q) 
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by multiplying Eq. (fT6l) by Q and integrating over dci and summing for i = a, b. 
It then follows that the term containing G vanishes on account of the subsidiary 
conditions, Eqs. (ICll) . One then seeks a solution satisfying the extremal condition 

consistent with Eqs. (ICip . after (pf^ is substituted by Eq. (fT2il . Using then the in- 
variance of the collision kernels under the exchange of velocities Cj ^ cfj, Ci ^ c'i and 
microscopic reversibility, one gets an equation for Q consistent with Eqs. (ICll ) which 
is of the form 

.2 5^ 



(0) = clO + ft^ [C (clO) + C (clO + c,0)] . 

Tfli 



\2kT 2 
If we now choose for (, 



M 



m=0 



and substitute back into the equation for -D(O), the left side is readily integrated 
leading to the result that 



5kTY: 



m 



■i (1) 
a) 



-SiBkTy^ \2( 



{K, L} 



where {K, L} is given by Eq. (IB14I1 . The remaining of the procedure is the same as 
in the previous case so that after maximizing the resulting set of algebraic equations, 
solving them, and separating real and imaginary parts [see Eqs. f lB3l) -( lB4ll ] one gets 
that 



Re [ai°)] = — (56.3 - 662x2 - 2.25x^) 



Im [4°)] = ^ (647x + 2.2x^) 



Re [a^^)] = ^ (37.5 + 2U7x'^ + 7.3x^) 



(B19) 
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Im [al^^ = (206a; + 26A9x^ + 9x^) 



Re 



r 



MiAi 



Im 



.(1) 



(1.12 + 121.2x2 + 154.4x^) 



(0.06x + 7x^ + 9x^) 



MiAi 



where Ai is defined in Eq. (El) and x = = 1.76 x lO^V^ if B is given in teslas. 
By a completely identical procedure, 



Re 



= (22.8 + 55.3x2 ^ o.l9x^) 

: ^ (234x + 331. + 1.12x^) 
Ai 

Re [4^)] = ^ (5.6 - 66.2x2 _ o.22x^) 
Im [4^)] = ^ (64.7x + 0.22x3) 



Im 



(B20) 



Re 



(1) 



Im 



r 



(1) 



(0.03-4.1x2 -0.22x^) 
: (.66x - 3.6x3) 



Notice that consistently in both cases, if i? = the real parts of these equations 
agree with the results given in Eqs. (IB19I1 and (IB20I1 . These equations are those 
needed to evaluate all the transport coefficients appearing in the text. 

* The authors are deeply grateful to Prof. Leonardo Dagdug for checking all the 
values of the collision integrals listed in this appendix and for his aid in all the rather 
involved algebraic steps involved in these calculations. 



Appendix C: SUBSIDIARY CONDITIONS 



In Sect. II it is asserted that Eq. (JH]) is a solution to the homogeneous part 
of Eq. ([7]). This statement has a deeper significance. In fact, as it is well known 
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any linear combination of the conserved variables nii, rriiVi and \mivf {i = a, h) 



is a solution to homogeneous part. This leads to a local distribution which has 
five undetermined parameters, all depending of position r and time t. To uniquely 
determine them one assumes that the five local variables n (r , t) , u (r , t) and e (r , t) 
are determined by ff'^ and not by the full distribution function fi {r,Vi,t). This is 
the way the local equilibrium assumption is forced into the kinetic formalism thus 
implying that, for a multicomponent system. 



where Lpl"' is the n*'^ coefficient in the Knudsen expansion of fi [Eq. ([9])]. Equations 
(IClll are the so-called subsidiary conditions that all ^^j's must obey. 




(CI) 
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Figure Captions 

Fig. 1. The Soret coefficients (5'||, solid; S±, dotted; Sg, dashed) for n — lO^^m"^ 
and T = lO^X. 

Fig. 2. The behavior of the three thermal conductivities {k\\, solid; k±, dotted; 
Ks, dashed) for n = lO^^m-^ and T = WK. 

Fig. 3. The behavior of the thermal diffusion coefficients D (D\\, solid; D^, dot- 
ted; Ds, dashed) consistent with the condition x < 1 for n = lO^^m"^ 
and T = WK. 

Fig. 4. The effective thermal conductivity k.*j_ compared with k.± for n = lO^^m"^ 
and T = WK. 

Fig. 5. The effective thermal conductivity compared with Kg for n — lO^^m'^ 
and T = WK. 
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